arXiv: 1502.01721 v2 [hep-ph] 2Jun2015 


Spot the stop with a 6-tag 


CERN-PH-TH-2015-019 

DFPD-2015/TH/03 


Gabriele Ferretti 

Department of Fundamental Physics, Chalmers University of Teehnology, 4 IS 96 Goteborg, Sweden 

Roberto Franceschini 

Theory Division, Physics Department, CERN, CH-1211 Geneva 23, Switzerland 

Christoffer Petersson 

Department of Fundamental Physics, Chalmers University of Technology, 4^2 96 Goteborg, Sweden 
Physique Theorique et Mathematique, Universite Libre de Bruxelles, C.P. 231, 1050 Brussels, Belgium and 

International Solvay Institutes, Brussels, Belgium 

Riccardo Torre 

Dipartimento di Fisica e Astronomia, Universita di Padova, and INFN Sezione di Padova, Italy 

The LHC searches for light compressed stop squarks have resulted in considerable bounds in the 
case where the stop decays to a neutralino and a charm quark. However, in the case where the 
stop decays to a neutralino, a bottom quark, and two fermions via an off-shell IF-boson, there 
is currently a significant unconstrained region in the stop-neutralino mass plane, still allowing for 
stop masses in the range 90-140 GeV. In this paper we propose a new monojet-like search for light 
stops, optimized for the four-body decay mode, in which at least one 6-tagged jet is required. We 
show that, already by using the existing 8TeV LHC data set, such a search would cover the entire 
unconstrained region. Moreover, in the process of validating our tools against an ATLAS monojet 
search, we show that the existing limit can be extended to exclude also stop masses below 100 GeV. 
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Introduction .—The top quark gives rise to the lead¬ 
ing quantum correction that destabilizes the electroweak 
scale in the Standard Model (SM). One way to solve 
this so-called hierarchy problem is to extend the ordi¬ 
nary spacetime symmetries by supersymmetry (SUSY) 
and introduce new physics at a low scale in the form of 
superpartners of the SM particles. In order to cancel the 
leading quantum correction, the superpartner of the top 
quark, the stop, should have a mass of the order of the 
electroweak scale and hence be observable at the Large 
Hadron Collider (LHC). Light stops have been subject 
to intense recent studies both in the theory [1-36] and 
experimental [37-50] communities. 

Taking a simplified model approach, a key strategy to 
test R-parity conserving SUSY is to consider only the 
lightest stop mass eigenstate U, decaying to the light¬ 
est superpartner (LSP), the neutralino x®, and taking 
all other superpartners to be sufficiently heavy and ef¬ 
fectively decoupled. In such a simplified model, where 
the stop mass and the neutralino mass are the 
only two parameters, the only relevant SUSY production 
mode is stop pair production, for which the cross-section 
is determined by . 

In the case where the mass difference between the 
stop and the neutralino is larger than the top mass, 
— m^o > m*, each of the pair produced stops 
decays 2-body via an on-shell top, U —f t Xi- In this case, 
for neutralino masses below around 250 GeV, the cur¬ 
rent LHC limits exclude stop masses below 600-750 GeV 


[41, 50], while for larger neutralino masses there are no 
bounds. 

For mass splittings in the range mw + mt < Am < mt , 
the stop decays 3-body via an off-shell top, ii^bW x°, 
and the stop mass limits reach up to around 
200-300 GeV [38, 41, 44]. It should be noted that the 
limits in the different mass splitting regions are not con¬ 
tinuously connected to each other and close to the mass 
thresholds at Am ^ m* and Am ~ m\Y + mt, the bounds 
become weak or disappear. 

For even smaller mass splittings, Am K-mw+iTib, the 
stop can have two different decay modes, either the 
4-body decay, via an off-shell W, ii^bff'xi, or 
the 2-body decay U—^-cxi- While the former decay 
mode is phase space suppressed, the latter is 1-loop 
and Cabibbo-Kobayashi-Maskawa (CKM) suppressed, 
see e.g. Refs. [21, 34, 51, 52] for discussions. The branch¬ 
ing ratios for these two competing decay modes depend 
strongly on the flavor structure of the squark soft masses, 
such as the off-diagonal stop-scharm mixing mass term, 
as well as the masses of the superpartners that enter the 
loop. In order to be able to restrict to the minimal set of 
parameters, it is customary to consider these two decay 
modes separately, and assume 100% branching ratio in 
each case. 

In the small mass splitting case, several LHC searches 
have been optimized for the 2-body charm decay mode. 
Under the assumption BR(ti —>• cxi)=l, stop masses be¬ 
low 250-300 GeV have been excluded [42, 45]. In con- 
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trast, no existing LHC analysis has been optimized for 
the 4-body stop decay mode. Currently, even though 
the 4-body case is partly covered by the ATLAS searches 
[41, 42], there exists a significant unconstrained region 
in the stop-neutralino mass plane, still allowing for stop 
masses in the range 90-140 GeV (for neutralino masses 
below 60GeV). Since very light stops could be hiding 
there, it is important to find means to probe it. In this 
paper we show that, by augmenting an existing mono¬ 
jet search with a 6-tag requirement, this unconstrained 
region could be covered completely already with the ex¬ 
isting 8TeV LHC data set. 

Existing stop 4-body searches .—We start by reviewing 
the existing searches relevant for stops in the mass range 
nib < Am < for which BR(fi -^bff x5)=l. In 

contrast to the case where the stop decays via the 2-body 
charm decay mode, which has been probed by Tevatron 
[53], CMS [45] and ATLAS [42], neither Tevatron nor 
CMS have performed any search that places a bound in 
the 4-body decay case.^ Therefore we will focus our dis¬ 
cussion on two searches performed by ATLAS. 

The first search that places a limit in the stop 4-body 
decay case is a monojet search in which events are re¬ 
quired to contain at least one hard jet and a large amount 
of missing transverse energy [42]. The exclusion 

curve arising from this ATLAS search is indicated by the 
blue dashed curve in Figure 1 (left). As can be seen in the 
Figure, this search is most sensitive to the case where 
the mass splitting Am is small^ and the stop 4-body 
decay products are soft. The required hard jet arises 
from initial state radiation (ISR), against which the pair- 
produced stops recoil. The ISR jet boosts the two stops, 
which are no longer produced back-to-back, thereby in¬ 
creasing the Ipp in the event. 

The abrupt end of the blue curve at m^^ =100 GeV in 
Figure 1 (left) is simply due to the fact that ATLAS does 
not provide the limit for smaller stop masses. Given the 
requirements in this search, one would expect that the 
exclusion curve should continuously extend diagonally 
down to the left, reaching the LEP limit [54], which is 
indicated by the black dashed curve. In the Section Re¬ 
sults below we discuss this issue further and provide the 
expected limits for stop masses below 100 GeV. 

The second ATLAS search that places a bound in the 
stop case under consideration is a search in the final state 


^ Searches for the stop decay ii might have some sensitivity 

to the 4-body decay mode that we study. However, since the 
results are not presented in terms of the 4-body stop decay, we 
do not include them in our summary of the existing bounds. 

^ For Am smaller than about 20 GeV, the partial width for the stop 
4-body decay decreases to the point where the stop either decays 
via a displaced vertex or, if other decay channels are present, the 
4-body branching ratio is strongly suppressed. However, in the 
spirit of simplified models, we follow the same strategy of ATLAS 
and present our results for 100% 4-body decay branching ratio, 
assuming the stop to always decay promptly. 


with one lepton, jets and fp [41]. Since this search relies 
on the presence of a lepton, arising from the stop 4-body 
decay, it is most sensitive to the case where the mass 
splitting Am is at least sufficiently large to allow for a 
reconstruction of the lepton. The exclusion curve aris¬ 
ing from this search is indicated by the orange curve in 
Figure 1 (left). Note that this curve ends at neutralino 
and stop masses around 60 GeV and 110 GeV, respec¬ 
tively. Unlike the monojet search, in which the abrupt 
ending of the exclusion curve was due to the lack of sig¬ 
nal samples, the reason for this ending has a physical 
origin. As one moves diagonally down to the left, i.e. 
keeping fixed Am, the Ijdp spectrum becomes softer and 
the search looses sensitivity [55]. Hence, while the cross- 
section increases, as a consequence of the decreasing stop 
mass, the acceptance times efficiency decreases faster, as 
a consequence of the decreasing neutralino mass. 

Figure 1 (left) summarizes the current experimental 
status concerning searches for light stops that dominantly 
decay in a 4-body final state. We see that there is a tri¬ 
angular shaped unconstrained region for stop masses in 
the range 90—140 GeV, with boundaries given by the ex¬ 
clusion curves from LEP [54], the ATLAS monojet search 
[42] and the ATLAS 1-lepton search [41]. 

Proposed stop 4-body search .—The search we propose 
is a simple extension of the ATLAS inclusive monojet 
search [42] where, in addition, we require the presence 
of at least one 6-tagged jet. This 6-tag requirement is 
motivated by how efficiently it reduces the leading back¬ 
ground processes in Ref. [42], with respect to the stop 
4-body decay signal process, thus enhancing the sensi¬ 
tivity to light compressed stops. 

In order to display the potential gain in sensitivity that 
can be achieved, we consider the addition of a 6-tag to the 
definition of the signal region “Ml” of the ATLAS mono¬ 
jet search [42]. The Ml selections require the presence 
of at most three jets (including 6-jets) with pp > 30 GeV 
and \r]\ < 2.8. The leading jet pp must be greater than 
280 GeV and the Ifip in the event must be above 220 GeV. 
There is a veto on muons with pp > 10 GeV and \r]\ < 2.4 
and on electrons with pp > 20 GeV and \r]\ < 2.47. More¬ 
over, there is a condition requiring that the azimuthal 
angle Acj) between the Ipp-Mector and each of the jet 
PT-vectors should be larger than 0.4. The signal region 
that selects events according to these cuts is denoted by 
Ml, in accordance with the ATLAS notation. 

We remark that the selections that define the signal 
region Ml might not be optimal over the entire range of 
mass spectra that we consider. For instance, one could 
apply the same argument of adding 6-tagging information 
to the other signal regions defined in Ref. [42], in which 
harder Ifip cuts are employed. Moreover, as the mass 
splitting between stop and neutralino varies, one could 
optimize further the number of jets that defines the sig¬ 
nal region. This is done for instance in the Dark Matter 
search in Ref. [56] where two signal regions with at least 
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Background 

tt 

Z{-^uv) 


Dibosons 

Others 

Total 

Ml (ATLAS [42]) 
MU-fe-tag 

780 ± 73 

307 ± 57 

17400 ± 720 

261 ± 22 

14100 ± 337 

144 ± 7 

650 ± 99 

55± 17 

565 ± 301 

33450 ± 960 

767 ± 64 


TABLE I. Estimated numbers of background events with 20.3 fb ^ of 8TeV LHC data. 


one 6-jet and large Ipj, are employed. However, we find 
the Ml selections sufficient to make the case for the ad¬ 
dition of a 6-tag that we propose and hence, in this work, 
we choose to concentrate on this signal region. The ben¬ 
efit of Ml is that it is the signal region in Ref. [42] with 
the lowest cut and therefore it is the most promis¬ 
ing for low mass stops. Furthermore this signal region is 
defined more inclusively on the number of jets than the 
signal regions of Ref. [56] , which allows us to make more 
reliable predictions for the signal and background, and 
also to better reproduce the ATLAS results. 

The background and signal simulations were performed 
using MadGraphS [57], Pythia6 [58], FastJetS [59, 
60] and DelphesS [61] with the ATLAS standard detec¬ 
tor specification. Jets (including 6-jets) are reconstructed 
using the anti-A:* clustering algorithm [62] with a jet ra¬ 
dius parameter of 0.4. We have used the CTEQ6L1 PDF 
sets [63] and MLM jet matching [64, 65] throughout the 
analysis. 

The leading SM backgrounds for our analysis are the 
tt, Z{-^vv)+ieis, W{-^^v)+iets (where £ = e,^,T) and 
diboson processes. For the normalization of the cross 
sections of these processes at 8TeV LHC, we have used 
the same theoretical predictions as ATLAS in Ref. [42], 
obtained from Refs. [66-72]. We have validated our simu¬ 
lations against the Ml signal region of the ATLAS mono¬ 
jet search. Comparing our estimates for the background 
rates with the corresponding ATLAS estimates in Table 
VHI of Ref. [42] , we find that the central values from our 
simulations are within 20% of the corresponding ATLAS 
numbers. This agreement is quite satisfactory and makes 
us confident that reliable conclusions can be drawn using 
our tools. To further improve our estimates, we normal¬ 
ize our predictions for all the leading background pro¬ 
cesses to exactly match ATLAS background estimates, 
which, for the Z-|-jets and W-l-jets samples, have the ex¬ 
tra advantage of incorporating data-driven re-weighting 
of the simulated events. The expected backgrounds are 
summarized in Table I written as R ± SB, where B is 
the central value and SB the la error. For all the back¬ 
grounds shown in Table 1 we quote Icr errors from AT¬ 
LAS. 

The crucial ingredient in our analysis is that we ex¬ 
tend the Ml set of selections by adding the requirement 
that the events passing the Ml cuts must contain at least 
one 6-tagged jet with px in the range 30—300GeV, and 
with |77|<2.5. The 6-tagging is parametrized through 
DelphesS [61] using a “mild” working point character¬ 
ized by a light flavor jet rejection of 1/1000 and an effi¬ 


ciency for actual 6-quarks of around 0.4 for central jets. 
We emphasize that for the px range that we have cho¬ 
sen for the 6-tagged jet, the calibration of the 6-tagging 
algorithms in the ATLAS experiment is data-driven, min¬ 
imizing the systematic uncertainties on 6-tagging coming 
from Monte Carlo simulations [73] . We denote this signal 
region by “Ml-|-6-tag”. 

In the row Ml-|-6-tag in Table 1 we show the number 
of background events we get in the proposed signal re¬ 
gion. The comparison of the events expected for the Ml 
and for the Ml-|-6-tag signal regions gives an idea of the 
power of the 6-tag requirement to enhance the sensitiv¬ 
ity to the signal. In fact, the stop signal is expected to 
behave similarly to the ti background, hence to be only 
mildly reduced by the 6-tag requirement. Typical sig¬ 
nal efficiencies to the 6-tag requirement are around 20%. 
In contrast, we observe a great suppression of the W 
and Z boson backgrounds, which are the leading back¬ 
grounds in the MI signal region. We stress that in our 
work the efficiency of signal and backgrounds to the 6- 
tag requirement is given by the parametrization imple¬ 
mented in Delphes3 [61], which is expected to be reli¬ 
able. In order to draw conclusions on safe grounds, we 
take the relative error of the Ml-|-6-tag background pre¬ 
diction to be twice the relative error quoted by ATLAS 
for the Ml region. This gives background relative errors 
slightly smaller than those quoted by ATLAS [42] in the 
signal regions involving c-tags. Given that 6-tags are un¬ 
der better experimental control than c-tags, we expect 
our estimate of the uncertainties to be fair. 

The signal process has been simulated for a grid of 
points with varying from 70 GeV to 250 GeV, and 
m^o varying from 0 to 200 GeV, in steps of 10 GeV. The 
only relevant SUSY production mode is stop pair pro¬ 
duction, with cross section given in Table H. For stop 
masses in the range 100-250 GeV we have used the next- 
to-leading-order (NLO) plus next-to-leading-logarithm 
(NLL) stop cross sections used by ATLAS and given 
by the LHC SUSY Cross Section Working Group [74], 
to which an uncertainty of around 16% is assigned. For 
the stop mass points below 100 GeV, which are not given 
by ATLAS, we computed the NLO cross sections using 
Prospino [75] and normalized them with the available 
ATLAS cross sections for stop masses above 100 GeV. We 
compared the expected number of signal events we ob¬ 
tained to the corresponding numbers reported by ATLAS 
[42] in several points and we found a systematic overshoot 
of around 20%. Therefore we normalized the expected 
number of signal events by decreasing it of 20% to match 
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FIG. 1. Left panel: Summary of the existing limits in the stop-neutralino mass plane, as well as our exclusion curve for the 
Ml signal region of the ATLAS search [42], in which we also extend the range of excluded stop masses below 100 GeV. Right 
panel: Exclusion arising from our proposed search Ml+b-tag, indicated by the red solid curve. The red dashed curve denotes 
the change caused by a 20% increase of the total background error. 


the ATLAS central values. 

The efficiencies in the analysis are rather small, re¬ 
quiring the generation of a large number of Monte Carlo 
events. For the backgrounds in the Ml case we are able to 
generate a sufficient amount of fully jet-matched events 
to keep the statistical error below 10%. However, for the 
signal, given the amount of points in the grid we aim to 
cover, we do not have the computer resources to generate 
fully matched events at a similar level of statistical uncer¬ 
tainty. We solve this problem by performing the analysis 
in two steps. Since all events that pass the cuts of our 
signal regions contain a hard jet, for each point of the 
grid, we generated both the exclusive zero-jet leading- 
order (LO) stop pair production process pp^titi and 
the one-jet process pp^iiiij, with jet pT> 200 GeV. 
The ratio of these two cross sections is used to calculate 
the efficiency of the px >200 GeV cut. The unmatched 
one-jet sample was then used for the grid in the analy¬ 
sis. We have checked the validity of this procedure in a 
few points by generating fully matched inclusive samples, 
with sufficient statistics, and we found the two proce¬ 
dures to be in agreement within one statistical standard 
deviation. 

Results .—In order to estimate whether a point in the 
stop-neutralino mass plane is excluded, we compute the 
number of expected events for that point in a given sig¬ 
nal region. We expect that the experiments can put 
a 95% C.L. exclusion for those mass points that yield 
a number of signal events greater than V 95 = 1.96 6B, 
where SB is the total error in Table I. 

Following this limit-setting procedure we start by cal¬ 
culating the exclusion we obtain using the Ml signal re¬ 


gion. The resulting exclusion curve is given by the red 
solid curve in the left panel of Figure 1. We remark that 
our exclusion curve follows quite closely the one given by 
ATLAS, which further validates our procedure. It should 
be noted that our curve extends down to the LEP bound 
for m^a > — 40 GeV, thereby covering part of the un¬ 

constrained region between the blue ATLAS curve and 
the black LEP curve. Hence, by considering stop masses 
below 100 GeV, the existing ATLAS bound arising from 
the Ml signal region can be extended. 

Following the same limit-setting procedure for the sig¬ 
nal region Ml-|- 6 -tag, we obtain the red solid exclusion 
curve in the right panel of Figure 1. This is the main 
result of this paper. The dashed red curve corresponds 
to a 20 % increase of the total error on the background. 
The comparison of the two red curves gives an idea of the 
sensitivity of our result to a) the uncertainty associated 
with the signal cross sections, h) possible contributions 
from sub-leading backgrounds not evaluated for Ml-l- 6 - 
tag. We see that our proposed search Ml-|- 6 -tag cov¬ 
ers the entire unconstrained region. Moreover, it slightly 
extends the existing LHC limits for stop mass around 
200 GeV. 

It is worth mentioning that our simulations suggest 
that further sensitivity is gained by removing the 
condition, but keeping the same 6 -tag requirement as in 
Ml-|- 6 -tag. In the Ml signal region, the condition is 
introduced to reduce the pure-QCD multijet background, 
for which the fj. originates from jet mismeasurements. 
However, the 6 -tag requirement can be seen as an al¬ 
ternative to the A(/) cut since it is expected, already by 
itself, to dramatically reduce the multijet background. 
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mi [GeV] 

70 

80 

90 

100 

110 

120 

130 

140 

150 

160 

170 

180 

190 

200 

210 

220 

230 

240 

250 

a (it) [pb] 

2797 

1550 

912 

560 

362 

240 

163 

113 

80.3 

58.0 

42.6 

31.9 

24.2 

18.5 

14.3 

11.2 

8.78 

6.97 

5.58 


TABLE II. Inclusive stop pair production cross-sections used in the analysis. 


With our simulation tools, the estimation of the multi¬ 
jet background would not be reliable, therefore we do 
not attempt to estimate the gain in sensitivity to light 
stops that could be achieved by employing a looser A(^ 
cut. Instead we content ourselves with simply encour¬ 
aging the experimental collaborations to also consider a 
signal region with a looser A(j) cut. 
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